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Abstract 



In this letter we give the equation of state of QCD at finite temperatures and densities. The recently proposed overlap 
improving multi-parameter reweighting technique is used to determine observables at nonvanishing chemical potentials. 
Our results are obtained by studying nj =2+1 dynamical staggered quarks with semi-realistic masses on Nt — 4 lattices. 

wo: 

^ ' According to the standard picture of QCD, at high temperatures and/or high densities there is a change from a state 
ominated by hadrons to a state dominated by partons. This transition happened in the early Universe (at essentially 
^Vanishing density) and probably happens in heavy ion collisions (at moderate but non- vanishing density) and in neutron stars 
C^at large density, for which a rich phase structure is conjectured [H @, §])• There are well established nonperturbative 
lattice techniques to study this transition at vanishing density, at which the equation of state were determined as a function 
' ^if i the temperature. Due to the sign problem (oscillating signs lead to cancellation in results, a phenomena which appears in 
O^iany fields of physics) nothing could have been said about the experimentally/observationally relevant case at non- vanishing 
(---densities. By using our recently proposed lattice technique || we calculate the equation of state at non- vanishing temperatures 
(^ind densities. In this letter we present the technique and show the results. A detailed version of the present work will be 
OQmblished elsewhere ||. 

' QCD at nonzero density can be studied by introducing the chemical potential (/x). Such a theory is easily formulated on 
("-^pace-time lattice @, However; standard, importance sampling based Monte-Carlo techniques can not be used at fi ^0. 

**T*Up to now, no technique was suggested capable of giving the equation of state (EOS) at non-vanishing fi, which would be 
essential to describe the quark gluon plasma (QGP) formation at heavy ion collider experiments. Results are only available 
ifor /z=0 (Til 11, O) at non- vanishing temperatures (T) (a recent review is Ref. pT 



Oh, The overlap improving multi-parameter reweighting [^| opened the possibility to study lattice QCD at nonzero T and //. 

irst one produces an ensemble of QCD configurations at /x=0 and at T^O. Then the Ferrenberg-Swendsen type reweighting 
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ctors 14 1 of these configurations are determined at \i ^ and at a lowered T. The idea can be easily expressed in terms of 
^he partition function 



Z{ji,p) = J T>Uexp[-Sg{P,U)}detM(Li,m,U) = 
J VU exp[-S g (/3 ,U)]detM(fi = 0,m,U) 

where S g is the action of the gluonic field (U), (3 = 6/g 2 fixes the coupling of the strong interactions (<?). Note that for a 
given lattice T is an increasing function of (3. The quark mass parameter is in and det M comes from the integration over the 
quark fields. At nonzero \x one gets a complex det M which has no probability interpretation, thus it spoils any importance 
sampling. Therefore, the first expression of eq. ([!]), U ^0, is rewritten in a way that the second line of eq. (0) is used as an 
integration measure (at /i=0, for which importance sampling works) and the remaining part in the curly bracket is measured 
on each independent configuration and interpreted as a weight factor {w(/3, (A,m, U)}. In order to maximise the accuracy of 
Z the reweighting is performed along the best weight lines on the \i-(3 plane (or equivalently on the fi-T plane). These best 
weight lines are defined by minimising the spread of log w. 

Similar reweighting can be done in the mass parameter, too. Using the above technique, transition (or hadronic/QGP) 
configurations are reweighted to transition (or hadronic/QGP) configurations as illustrated by Fig. Q thus, a much better 
overlap can be obtained than by reweighting pure hadronic configurations to transition ones as done by single /x-reweighting 
|15|| . As we emphasised, the technique works for temperatures at, below and above the transition temperature (T c ). By using 
the reweighting technique, the phase diagram Q and the location of the critical endpoint [ fl6| was given (for other approaches 
see e.g. [0 Using a Taylor expansion around /x=0, T^O for small jj, is a variant of our multi-parameter reweighting 
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Figure 1: The best weight lines on the plane, along which reweighting is performed. In the middle we indicate the 
transition line. Its first dotted part is the crossover region. The blob represents the critical endpoint, after which the transition 
is of first order. Below the transition line the system is in the hadronic phase, above the transition line we find the QGP. 
The integration paths used to calculate the pressure are shown by the arrows along the /i = axis and the best weight lines. 




Figure 2: The pressure normalised by T 4 as a function of T/T c at /i = (to help the continuum interpretation the raw 
lattice result is multiplied with c p ). The continuum SB limit is also shown. 

method, which can be used to determine thermal properties ]19[ |. A completely different method, analytic continuation from 
imaginary fi, confirmed the result of the reweighting technique on the phase diagram pt| . 

In the present analysis we use N t ■ finite T lattices with 7V t =4 and iV s .=8,10,12 for reweighting and we extrapolate to 
the thermodynamical limit using the available volumes (V). At T=0 lattices of 24 ■ 14 3 are taken for vacuum subtraction and 
to connect lattice parameters to physical quantities. 16 different f3 values are used, which correspond to T/T c = 0.8, . . . , 3 
and at T c to a lattice spacing (a) of «0.28 fm. We use 2+1 flavours of dynamical staggered quarks. Only fj, values of the light 
quarks are studied. While varying [3 (thus the temperature) we keep the physical quark masses constant at m u d ~ 65 MeV 
and m s ~ 135 MeV (the pion to rho mass ratio is to w /to p «0.66). From now on we usually omit the different quark mass 
indices. 

The determination of the equation of state at fj, ^0 needs several observables O at non-vanishing (i values. This can be 
calculated by using the weights of eq. (|l|) 

-p* lR \ EMA Mi rn, U)}Q(P, fM, m, U) 

£>(/?, /i, m) = — . 2 

2^{w{p,fi,m,U)} 

We use the following notation for subtracting the vacuum term: (0(/3, /i, Tn))= 0(/3, [a, m) T _ i() — 0((3, fi — 0, m) T _ . 

The pressure (p) can be obtained from the partition function as p=T ■ d log Z/dV which can be written as p—(T/V) ■ log Z 
for large homogeneous systems. On the lattice we can only determine the derivatives of log Z with respect to the parameters 
of the action (f3, m, n), so p can be written as an integral[|T): 

P 1 Lur \ (/ d{\ogZ) \ / d{\ogZ) \ / d{\ogZ) \\ 
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Figure 3: Ap = p(/z ^ 0, T) - p(/x = 0,T) normalised by T 4 as a function of T/T c for /i B =100, 210,330,410 MeV and 
Hb =530 MeV (from top to bottom; to help the continuum interpretation the raw lattice result is multiplied with c M ). 
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Figure 4: Ap of the interacting QCD plasma normalised by Ap of the free gas (SB) as a function of T/T c for 
^s = 100,210,330,410 MeV and p B =530 MeV. The result is essentially /^-independent. 

The integral is by definition independent of the integration path. The chosen integration paths are shown on Fig [l]. 

The energy density can be written as £ = (T 2 /V) ■ d(\og Z)/dT + (/iT/V) ■ d(log Z)/d/j,. By changing the lattice spacing 
T and V are simultaneously varied. The special combination e — 3p contains only derivatives with respect to a and /i: 

e — 3p a d\og(Z) 

T 4 ~ ~ da 

The quark number density is n — (T/V) ■ d log(Z) /dfi which can be measured either directly or obtained from the pressure 
(note, that the baryon density is ns=n/3 and the baryonic chemical potential is /zs=3/i). 

We present lattice results on p(/x = 0,T), Ap(p,T) = p(/x ^ 0,T) - p(/i = 0,T), e(fi, T)-3p(/x, T) and n B {p,T). Our 
statistical errorbars are also shown. They are rather small, in many cases they are even smaller than the thickness of the lines. 
On the figures we multiply the lattice results with the dominant correction factors between 7V t =4 and the continuum in the T — > 
oo case: c p — p{ji = 0, T — > oo, continuum) /p(/i = 0, T — > oo, N t = 4)=0.518 and c p = Ap(/i, T — > oo, continuum) /Ap(/i, T — > 
oo,N t = 4)=0.446. The well-known continuum expressions in the T — > oo Stefan-Boltzmann (SB) case are p(fi = 0,T — > 
oo, continuum) = (16 + 21nf/2)7r 2 T 4 /90 and Ap(/i, T — > oo, continuum) = nf/i 2 T 2 /2 + 0(/j, 4 ). This way the results presented 
on the figures might be interpreted as continuum estimates and could be directly used in phenomenological applications. 
Direct Nt=4 lattice results can be obtained by dividing the presented values with these c p ,c M correction factors. 

Fig. H shows the pressure at /i=0. On Fig. || we present Ap/T 4 for five different /i values. Fig. || gives Ap(/z,T/T c ) 
normalised by Ap SB — Ap(/i,T — > oo). Notice the interesting scaling behaviour. Ap/Ap SB depends only on T and it is 
practically independent of /i in the analysed region. Fig. || shows e-3p normalised by T 4 , which tends to zero for large T. 
Fig. H gives the baryonic density as a function of T/T c for different /i-s. As it can be seen the densities exceed the nuclear 
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Figure 5: (e - 3p)/T 4 at ^ B =0, 210,410 MeV and ^ s =530 MeV as a function of T/T c (from bottom to top; to help the 
continuum interpretation the raw lattice result is multiplied with c p ) 




Figure 6: The quark number density normalized T 3 as a function of T/T c for /j,b =100,210,330,410 MeV and /is=530 MeV. 
(to help the continuum interpretation the raw lattice result is multiplied with c^). As a reference value the line starting in 
the left upper corner indicates the nuclear density. 



density by up to an order of magnitude. 

As an important finding we mention that in the present analysis the applicability of our reweighting method, the maximal 
\x value scales with the volume as /x max • a ~ (Nt ■ N 3 )~ - 25 . If this behaviour persists, one could -in principle- approach the 
true continuum limit (note, that for Nt — > oo a ~ 1/Nt ~ (Nt ■ Ng)~ ' 25 , thus /z max is constant). 

In this paper we studied the thermo dynamical properties of QCD at nonzero \i. We used staggered QCD with 2+1 quarks 
on iV t =4 lattices. We determined for the first time the equation of state at nonzero temperature and chemical potentials. 
Future analyses should be performed at smaller lattice spacings and quark masses. A detailed version of the present work 
will be published elsewhere ||. 
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